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MODELING  FAST  MODULATION  EFFECTS  IN 
CESIUM  ATOMIC  CLOCKS 


J.  P.  Hurrell,  W.  A.  Johnson,  S.  K.  Karuza,  and  F.  J.  Voit 
Electronics  Technology  Center 
The  Aerospace  Corporation 
P.O.  Box  92957 
Los  Angeles,  CA  90009 


Abstract 

A  theory  for  fast  modulation  effects  in  cesium  frequency  standards  has  been  developed  and 
expressions  derived  for  frequency  offsets  produced  by  Rabi  and  cavity  pulling  with  square-wave  phase 
modulation.  These  are  compared  with  equivalent  expressions  for  sine-wave  phase  modulation.  Ratios 
of  the  amplitudes  of  the  frequency  offsets  in  the  two  modulation  schemes  are  calculated  for  a  mono¬ 
velocity  beam  and  shown  to  be  in  agreement  with  measurements  performed  on  a  commercial 
standard.  This  analysis  allows  a  more  complete  evaluation  and  comparison  of  different  modulation 
schemes  for  cesium  atomic  beam  standards. 

INTRODUCTION 

A  number  of  different  modulation  schemes  are  used  to  implement  the  servo  loop  in  cesium  atomic 
clocks.  They  include  sine-wave  phase  (frequency)  modulation,  square-wave  phase  modulation  and 
square-wave  frequency  modulation.  The  choice  of  modulation  scheme  is  based  on  a  number  of  consid¬ 
erations,  e.g.,  ease  of  implementation,  optimization  of  cesium  flux  usage  and  Allan  variance,  minimized 
frequency  offsets.  Analysis  of  servo  loop  gain  has  been  performed  for  all  these  schemes  and  expres¬ 
sions  exist  to  model  the  gain  dependence  upon  loop  parameters.  However,  the  analysis  of  the  effects  of 
modulation  scheme  on  accuracy  and  long  term  stability  are  less  complete.  Sine-wave  and  square-wave 
frequency  modulation  effects  have  been  considered!1!  within  the  context  of  slow  modulation!2!,  but  fast 
modulation  transients  have  not  been  modeled  because  of  the  lack  of  an  atomic  description  of  those 
atoms  that  experience  the  step  change  in  microwave  phase  or  frequency.  A  theory  for  square-wave 
phase  modulation  is  outlined  in  this  paper. 

In  the  analysis  of  servo  loop  gain,  it  can  be  assumed  that  the  cesium  atoms  experience  a  constant  mi¬ 
crowave  frequency  and  the  effects  of  phase  modulation  only  perturb  the  relative  phase  of  the  precessing 
atoms  and  the  microwave  field  in  the  second  arm  of  the  Ramsey  cavity!3!.  In  particular,  the  transient 
effects  in  square-wave  phase  modulation  can  be  safely  omitted,  and  the  contribution  neglected  from 
those  atoms  which  experience  the  step  change  in  phase  during  their  resonant  interaction  with  the 
microwave  field.  The  error  incurred  in  the  calculation  of  servo  gain  slope  is  only  of  order  umr  where 
t  is  the  interaction  time  and  um  is  the  modulation  frequency.  Similarly,  by  neglecting  the  frequency 
variation  of  the  microwave  field,  small  errors  of  order  Aur  are  incurred  in  the  frequency  modulation 
calculations.  represents  the  depth  of  frequency  modulation. 

Spurious  servo  signals  are  introduced  by  neighboring  transitions  (Rabi  pulling)  and  cavity  mistuning 
(cavity  pulling)  and  affect  the  performance  of  cesium  frequency  standards.  Differences  in  the  velocity 
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distributions  of  |  +m  >  and  |  -m  >  bave  the  potential  to  generate  a  frequency  offset  by  causing 
a  sloping  background  to  the  |  0  >— »|  0  >  transition.  Both  ^-transitions  (Am  =  0)  and  weak  x- 
transitions  (Am  =  ±1}  may  contribute  to  this  background.  Since  the  background  and  the  servo 
gain  slope  have  different  microwave  power  dependencies,  the  resulting  offset  provides  a  means  for 
transforming  power  changes  to  frequency  changes!4'.  Cavity  pulling  can  render  cesium  clocks  sensitive 
to  both  microwave  power  level  and  temperature  changes.  Both  offsets  are  mediated  by  exactly  those 
contributions  which  could  be  safely  neglected  above. 

Apart  from  transient  effects  in  frequency  modulation  schemes,  the  unwanted  signals  can  be  calculated 
by  allowing  the  atoms  to  experience  a  constant  microwave  frequency  during  resonance,  but  one  that 
changes  in  time  before  succeeding  atoms  arrive.  This  preserves  the  conditions  for  standard  microwave 
resonance.  The  transient  conditions,  on  the  other  hand,  require  a  different  description.  This  paper 
describes  an  approach  to  include  the  transient  contributions  to  both  Rabi  and  cavity  pulling.  The 
results  comparing  sine-wave  and  square-wave  phase  modulation  offsets  are  compared  with  experiments 
performed  on  a  commercial  frequency  standard. 

RABI  PULLING 

First  order  time-dependent  perturbation  theory  is  used  to  calculate  the  background  signal  resulting 
from  transitions  which  are  far  off-resonance.  This  approach  can  be  easily  applied  to  any  modulation 
scheme  because  it  allows  the  phase  modulation  to  be  addressed  by  Fourier  decomposition.  Writing 

<  F  =  4,  m  |  X  |  F  =  3,  m!  >=  2 hcn  cos [ut  +  nuimt  +  6n)  (1) 

n 

for  a  particular  |  m  >— *|  m'  >  transition  leads  to  a  general  first  order  transition  probability,  W,  for  a 
single  velocity  given  by 


f  sin(A  +  no;m/2)rsin(A  +  mu?m/2)r  1 
l  (A  +  nwm/2)(A  +  rnwm/2)  / 

x  jcos[(m  -  n)wm[t  -  T(1  +  D/L)\  +  0n-  9„]  +  cos[(m  -  n)u/m(t  -  TD/L)  +  6m-  0„]}.(2) 


W  —  ^  ^  CnCm 

n  m 


The  notation  follows  [2]  where  2A  represents  how  far  the  microwave  frequency  is  above  resonance,  r 
is  the  transit  time  through  one  arm  of  the  Ramsey  cavity,  T  is  the  transit  time  between  arms  distance 
L  apart,  and  D  is  the  detector  distance  from  the  second  arm.  Velocity  averaging  has  been  presumed 
in  order  to  average  to  zero  the  interference  term  between  interaction  regions  (there  are  no  Ramsey 
oscillations  in  the  wings  of  a  transition!).  The  sin  x/x  terms  arise  from  assuming  perfect  square-wave 
amplitude  variation  of  the  microwave  field  on  passing  through  a  cavity  arm!4!.  The  expression  can  be 
simplified  by  linearizing  these  terms  with  respect  to  frequency  over  the  small  n,m  values  of  interest. 
Then  it  can  be  applied  to  a  particular  modulation  scheme  to  predict  the  fundamental  Fourier  signal 
component  appropriate  to  narrow-band  amplification  and  detection. 


a)  Sine-wave  modulation:  u(t)  =  u  +  Awcoswmt  with  Cn  =  cJn( and  6n  =  0, 


W  -  ^A(e*s“*;Arj{cos->-[t  -  T(1  +  D/L)\  +  cos„4  -  TD/l]}.  (3) 

This  result  has  been  derived  previously  by  De  Marchi!4!. 
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b)  Square-wave  modulation:  cos[wt  +  0(t)j,  with  0(t)  alternating  between  0  and  0  at  frequency  wm 
and  for  which 

co  =  ccos0/2,  0o  =  0/2;  C„  =  (2c/nx)  sin  0/2,  0n  =  0/2  for  n  odd. 

Only  terms  with  n  =  0,  ±1  contribute  and 

A w  times  the  phase  factors  in  Eq.(3)  may  be  interpreted  as  the  instantaneous  frequency  in  the  two  arms 
of  the  Ramsey  cavity.  (2wm/x)  times  the  similar  phase  factors  in  Eq.(4)  are  the  Fourier  components 
of  delta  function  responses  from  atoms  residing  in  the  two  interaction  regions  at  the  time  when  the 
step  changes  in  phase  occur.  These  expressions  for  W  have  to  be  averaged  over  velocity  distributions 
and  summed  over  separate  |  m  >— »j  m'  >  transitions  to  complete  the  calculation. 

These  signals  are  nullified  by  a  servo  loop  error  signal,  with  the  resulting  offset  determined  by  the 
servo  gain  slope.  The  gain  slopes  are  determined  for  the  usual  mono-velocity  expressions)3)  which 
must  also  be  separately  averaged  over  the  velocity  distribution  for  the  |  0  >— »|  0  >  transition.  All  the 
averaging  procedures  have  been  suppressed  for  the  sake  of  clarity  in  the  following  expressions.  The 
procedures  have  been  discussed  elsewhere)1). 


(a;  -  wc).  = 


(w  -  wc),,  = 


{gWAl}  CQS(<;mr/2cosa;m[t  -  T(l/2+  D/L)\ 
T  sin2  2crJi(fi)cosum[t  -  T{  1/2  +  D/L)\ 


with  f)  =  2  sinwmT/2. 

coswmr/2cosa;m[t  -  T(l/2  -  D/L)} 


e  ‘  j’(laal)  sin2  2cr  sin umT /2  cos  wm[t  -  T(l/2  +  D/L)] 

In  the  limit  of  small  modulation  index,  Aw  —»  0  and  (w  -  we),  =  (w  -  wc)#9.  At  finite  values  for 
Aw,  (w  —  w),  >  (w  —  we),9.  Exact  values  for  this  inequality  will  depend  upon  the  velocity  averages 
and  the  specific  values  chosen  for  servo  parameters.  However,  the  magnitudes  of  the  offsets  will  remain 
very  similar. 

The  perturbation  expansion  developed  for  W  is  in  powers  of  (c/A).  The  ratio  of  second  order  to  first 
order  terms  is  (c/A)2.  Noting  that  c  «  2.5  kHz  and  A  ss  20kHz,  it  is  apparent  that  the  series  is 
rapidly  convergent. 


CAVITY  PULLING 

Cavity  pulling  is  caused  by  the  asymmetry  of  the  microwave  spectrum  about  its  unperturbed  frequency, 
induced  by  a  difference  between  that  frequency  and  the  cavity  resonant  frequency.  The  fact  that 
square-wave  phase  modulation  leaves  the  frequency  unchanged  except  during  the  step  changes  in 
phase  leads  to  a  natural  modeling  approach  in  which  the  atomic  resonance  is  treated  exactly  before 
and  after  the  step  change,  while  the  transition  itself  is  treated  by  perturbation  theory.  Only  those 
atoms  residing  in  the  separate  arms  of  the  Ramsey  cavity  at  the  time  of  the  step  changes  are  affected. 
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If  the  transition  occurs  infinitely  fast,  i.e.  in  a  time  short  compared  to  the  precession  period  of  the 
atoms  in  the  microwave  field,  the  sudden  approximation  can  be  invoked  to  predict  that  the  transient 
itself  will  not  induce  any  atomic  response.  Of  course,  the  step  change  in  phase  will  result  in  a  servo 
signal  from  those  atoms,  which  will  merely  be  added  to  the  usual  servo  signal  arising  from  the  greater 
number  of  atoms  in  flight  between  the  two  arms  of  the  cavity.  However,  the  finite  bandwidth  of  the 
microwave  cavity  prevents  infinitely  fast  changes  in  phase  from  occurring  and  causes  a  breakdown  of 
the  sudden  approximation.  Now  the  transient  frequency  changes  are  finite.  Moreover,  the  positive  and 
negative  frequency  transients  remain  symmetric  only  if  the  microwave  frequency  is  located  at  cavity 
resonance.  It  is  this  asymmetry  which  records  the  cavity  offset  and  induces  a  cesium  beam  signal  at 
the  modulation  frequency.  This  signal  is  derived  from  fast  responses  at  the  beginning  and  end  of  each 
half-cycle  of  the  servo  error  signal  and  may  be  written  as 


_  sjn2  2cr|(2cr)/tan(2cr)|  cos  wmT/2coswm[t  -  T(l/2  +  D/L)}  ■  (7) 

When  this  is  combined  as  before  with  the  servo  gain  slope,  the  offset  (a;  -  wc)  can  be  calculated.  For 
a  mono-velocity  beam,  it  reduces  to 


-8Q7Lum  cos  ujmT /2 
u%T  sinwmT/2 


The  similar  result  for  sine-wave  modulation!1!  may  be  written 


/u  —  u>c  \ 

\  W0  )  $ 


~4Q\Au  coswmT/2 


ji  + W)--W) 


(2er)/tan(2cr) 


(8) 

(9) 


In  the  limit  Aw  — »  0,  the  two  expressions  become  equal,  and  cavity  pulling  effects  are  the  same  for 
sine-wave  and  square-wave  modulation.  The  sine-wave  frequency  offset  is  reduced  as  Aw  increases 
from  zero. 


EXPERIMENTAL  RESULTS 

Sine-wave  and  square-wave  modulation  were  compared  by  measuring  the  frequency  of  an  atomic 
standard  for  both  modulation  schemes,  under  as  identical  conditions  as  possible.  This  allows  the 
differences  in  modulation  schemes  to  become  apparent  but  does  not  attempt  to  optimize  performance. 
The  square-wave  modulation  depth  9  was  held  close  to  90  degrees  and  the  modulator  drive  signal 
was  passed  through  a  narrow  band,  unity  gain  filter  and  phase  shifter  to  generate  the  sine-wave 
modulation.  This  produced  a  modulation  index  (Aw/wffl)  of  unity  and  allowed  phase  adjustment  to 
maximize  the  signal.  The  modulation  frequency  was  close  to  maximizing  the  gain  slope  of  the  servo 
loop  with  u>mT  fa  *  and  f)  ta  2. 

Rabi  pulling  was  separately  identified  by  measuring  the  frequency  as  a  function  of  C-field  at  two 
microwave  power  levels;  the  optimum  power  maximized  the  average  beam  current,  and  the  higher 
power  was  3  dB  above  optimum.  Figure  1  shows  the  resulting  fractional  change  in  frequency  with 
power  level  for  both  modulation  schemes.  Cavity  pulling  was  determined  by  measuring  the  frequency 
for  cavity  tuning  on  resonance  and  ±12  MHz  off-resonance.  In  the  data  reduction,  it  was  assumed 
that  the  Rabi  pulling  was  independent  of  cavity  tuning,  and  the  cavity  pulling  was  independent  of 
C-field.  Three  cavity  settings,  two  power  levels  and  two  modulation  schemes  lead  to  a  total  of  twelve 
curves  of  frequency  versus  C-field.  The  Rabi  pulling  could  be  reduced  to  a  single  curve  at  optimum 
power  by  the  following  relationships: 


Fig.  1.  C-field  dependence  of  the  fractional  change  in  frequency  caused  by 
a  3  dB  increase  in  microwave  power  above  optimum  (maximum  dc  signal) 
for  sine-wave  and  square-wave  phase  modulation. 


(u>  -  <jje)Mq  =  0.94 (w  —  ue)t  for  both  power  levels,  (10) 

and  (w  -  we)(P  +  3dB  =  1.85(w  -  ue)[P)  for  both  modulation  schemes.  (11) 


This  curve  is  shown  in  Fig.2  and  illustrates  that  Rabi  pulling  effects  are  small  in  this  standard.  Zero 
frequency  offset  may  be  arbitrarily  chosen  for  the  highest  C-field,  but  is  truely  zero  only  at  even  higher 
fields.  The  ratio  0.94  between  square-wave  and  sine-wave  amplitudes  is  closer  to  1  than  predicted 
from  Eqs.(5)  and  (6).  In  the  absence  of  velocity  averaging,  this  ratio  is  predicted  to  be  J\(2)  =  0.58. 
Also,  the  ratio  between  high  power  and  optimum  power  is  predicted  to  be  2.5,  slightly  higher  than 
the  measured  value  of  1.85.  It  is  reasonable  to  anticipate  that  velocity  averaging  will  improve  the 
agreement  between  theory  and  experiment.  The  approximately  linear  dependence  on  power  arises 
because  the  off-resonance  transition  probabilities  increase  linearly  with  power,  whereas  the  resonant 
transitions  are  already  optimized. 

Table  I  lists  the  measured  offset  frequencies  versus  cavity  tuning.  Systematic  errors  of  unknown  origin 
exceed  the  standard  deviation  (2  x  10~13)  in  the  frequency  measurements.  Zero  offset  was  arbitrarily 
chos ; <1  to  be  that  for  the  sine-wave  modulation  case  at  optimum  power  and  cavity  tuned  on  resonance. 
Eqs.(8)  and  (9)  would  predict  a  mono-velocity  offset  ratio 


(faJ  -  Ue)$q  _  2Ji(2) 

(w  -  ue),  1  +  2)  -  J2(2) 


(12) 
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TABLE  I. 


Fractional  frequency  offsets  (parts  in  1012)  determined  from  measurements 
at  two  microwave  power  levels  for  both  sine-wave  and  square-wave  modula¬ 
tion.  The  cavity  was  tuned  to  the  cesium  frequency  Jo  and  to  Jo±  12MHz. 
The  cavity  pulling  is  normalized  to  the  frequency  for  sine-wave  modu¬ 
lation  at  optimum  power  (maximum  dc  signal)  and  cavity  tuned  to  Jo- 
Standard  deviation  =  3  x  10“ 1S. 


Jo  +  12  MHz  F0  Jo  -  12  MHz 


optimum  power  sine-wave  -0.9  0  3.6 

square-wave  -3.6  -1.2  2.2 

optimum  +  3dB  sine-wave  -3.6  -0.2  8.5 

square-wave  -8.2  -1.0  6.9 


The  measured  ratio,  obtained  by  taking  the  difference  of  the  offsets  at  ±  12MHz,  is  1.3  at  both  power 
levels.  The  power  dependence  of  the  offsets  is  contained  in  the  usual  expression  (2cr)/tan(2cr)  which 
is  zero  for  4cr  =  x  and  rises  to  —1.7  as  the  power  is  raised  3  dB.  Clearly,  the  optimum  power  did  not 
produce  zero  cavity  pulling,  but  the  increase  by  a  factor  2.7  on  raising  the  power  3  dB  would  imply 
that  the  nominal  power  was  within  1.6  dB  of  that  necessary  for  zero  offset. 
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CONCLUSIONS 


Our  measurements  show  that  sine-wave  and  square-wave  phase  modulation  yield  similar  values  for 
Rabi  and  cavity  pulling  effects  in  cesium  frequency  standards.  The  result  is  consistent  with  the  theory 
presented  in  this  paper.  It  predicts  that,  in  the  limit  of  small  frequency  modulation  depth,  sine-wave 
modulation  produces  the  same  frequency  offsets  as  square-wave  modulation.  At  higher  modulation 
depth,  however,  Rabi  pulling  tends  to  be  larger  and  cavity  pulling  smaller  for  the  sine-wave  modulation 
scheme.  Exact  predictions  will  require  knowledge  of  specific  servo  parameters  and  averaging  over  beam 
velocities. 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs, 
specializing  in  advanced  military  space  systems.  The  Corporation’s  Technology  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staffs 
wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and  program 
support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are  provided  by 
these  individual  Technology  Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics,  VLSI 
reliability,  compound  semiconductors,  radiation  hardening,  data  storage  technologies, 
infrared  detector  devices  and  testing;  electro-optics,  quantum  electronics,  solid-state 
lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical  laser 
development,  optical  resonators,  beam  control,  atmospheric  propagation,  and  laser 
effects  and  countermeasures;  atomic  frequency  standards,  applied  laser  spectroscopy, 
laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent  imaging,  solar 
cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of  new 
materials:  metals,  alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of 
carbon;  development  and  analysis  of  thin  Aims  and  deposition  techniques; 
nondestructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  development  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch 
vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  and 
vulnerability  assessment;  contamination,  thermal  and  structural  control;  high 
temperature  thermomechanics,  gas  kinetics  and  radiation;  lubrication  and  surface 
phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing 
using  atmospheric  radiation;  soiar  physics,  infrared  astronomy,  infrared  signature 
analysis;  effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation;  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radiative 
signatures  of  missile  plumes,  and  sensor  out-of-field-of-view  rejection. 


